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1.  INTRODUCTION 


1.1  Manufacturing  Issues  and  Plv  Waviness.  The  quality  and  in-service  performance  of  composite 
structures  is  strongly  influenced  by  the  manufacturing  process.  From  autoclave  curing  of  flat  panel 
laminates  to  injection  molding,  the  manufacturing  process  wiU  inevitably  influence  the  perfomtance  of  the 
finished  component.  Quite  often,  "manufacturing  defects"  cannot  be  entirely  eliminated  and.  tl^rcfore. 
must  be  tolerated.  When  designing  composite  stmetures,  it  is  important  to  have  a  tlwrough  understanding 
of  the  influence  typical  manufaemring  defects  can  have  on  performance. 

A  frequently  encountered  manufacturing  defect  found  in  thick-section  composite  structures  is  ply 
waviness.  Ply  waviness  is  shown  schematically  in  Figure  1  and  represents  a  perturbation  in  ply  orientation 
from  the  ideal  case.  The  manufacture  of  composite  cylinders,  for  example,  is  generally  performed  with 
the  filament  winding  process  where  impregnated  fiber  bundles  or  tows  are  wound  onto  a  mandrel  of  the 
desired  shape.  Ideally,  fiber  placement  of  the  tows  is  such  that  they  confonn  to  the  mandrel  surface  in 
the  circumferential,  axial,  or  some  desired  off-axis  orientation.  Undulation  of  the  individual  plies  or  tows 
may  result  from  the  manufacturing  process.  Ply  waviness  can  significantly  influence  the  performance  of 
the  cylinder  by  reducing  its  stiffness  and  strength. 


Figure  1.  Plv  Waviness  in  a  Filament  Wound  Composite. 
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Currently,  the  manufacture  of  thennosetting  ccxnposite  cylinders  is  acctnni^shed  by  first  winding  the 
entire  cylinder  followed  by  a  post-curing  and  consolidation  pio<^.  This  process  results  Ui  {^y  waviness 
through  the  thickness  of  the  cylinder  in  both  axially  and  circumferentially  oriented  (4ies.  In  contrast, 
filament-wound,  thermoplastic,  thick-walled  cylinders  can  be  manufactured  with  an  "on-line"  ccmsolidation 
process,  reducing  ply  waviness  in  the  circumferential  direction. 

The  effects  of  ply  wavine^  in  composite  cylinders  (or  composite  structures  in  general)  can  be  studied 
by  considering  the  effects  of  ply  waviness  in  typical  sublaminate  (i.e.,  fiat)  configurations.  Such  an 
apinoach  has  been  adopted  in  this  work  where  a  generic  {90/0/90]  sublaminate  construction  with  a  wavy 
[0]  ply  is  investigated.  In  this  report,  an  analytic  model  is  developed  to  investigate  the  influence  of  ply 
waviness  on  the  stiffness  and  strength  of  composite  laminates.  The  model  predicts: 

•  elastic  properties  and  thermal  expansion  coefficierus  as  a  function  of  the  constituent  material 
properties  and  degree  of  undulation  (ply  waviness) 

•  individual  ply  stresses  for  prescribed  mechanical  and  thermal  load  cases 

•  the  strength  reduction  associated  with  ply  waviness. 

Parametric  studies  are  presented  which  quantify  the  synergistic  effects  of  ply  waviness  on  stiffness  and 
strength. 

The  theory  developed  in  this  work  represents  a  useful  design  tool  from  which  important  manufacturing 
and  design  trade-offs  can  be  made.  The  results  presented  indicate  that  the  degree  of  waviness  can 
significantly  reduce  the  stiffness  and  strength  of  the  structure. 

TI%  capabilities  developed  in  this  work  for  predicting  the  reduction  in  stiffness  and  strength  due  to 
ply  waviness  also  represent  a  valuable  tool  for  optimum  structural  design.  For  example,  the  reduced 
stiffness  properties  predicted  can  be  employed  as  input  into  a  suitable  global  stress  analysis,  resulting  in 
a  more  realistic  effective  response  prediction  of  the  structural  response  as  compared  to  a  defect  free" 
analysis.  The  resulting  stress  fields  predicted  in  the  structure  could  then  be  used  in  conjunction  with  the 
failure  analysis  developed  in  this  work  to  predict  the  reduction  in  strength  associated  with  ply  waviness. 
This  design  methodology  ultimately  results  in  a  more  realistic  prediction  of  strength  reduction  due  to  ply 
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waviness,  and  provides  valuable  insight,  in  a  consistent  manner,  into  the  dominant  failure  mechanisms 
characteristic  of  thick-section  composite  structures.  The  analytic  model  also  provictes  a  means  of 
quantifying  die  economic  trade-offs  that  exist  between  process  refit^ment  and  material  performance. 

1.2  Literature  Survey.  Tte  effects  of  the  waviness  of  the  ply  layers  on  the  mechanical  performance 
of  composite  strucnires  has  received  a  considerable  amount  of  attention  in  the  published  literature.  The 
focus  of  relevant  research  has  been  on  compression,  since  compression  performance  is  known  to  be  most 
influenced  by  waviness.  In  general,  waviness  in  composites  has  been  studied  on  two  different  levels: 
fiber  waviiress  on  tlK  fiber/inatrix  level  with  individual  fibers  embedded  in  a  surrounding  matrix  (i.e.,  in 
a  single  ply  or  lamina),  and  layer  waviness  of  a  complete  ply  or  plies  within  a  multidirectional  laminate. 
An  extensive  literature  review  on  waviness  in  composites  is  presented  by  Adams  and  Hyer  (1992).  A 
brief  suitunary  of  the  related  research,  however,  is  provided  below  for  completeness. 

The  following  literature  is  related  to  fiber  waviness  on  the  fiber/matrix  level.  Davis  (1975)  modeled 
fiber  waviness  in  boron/aluminum  composites  using  the  classic  microbuckling  model  (of  straight  fibers 
embedded  in  an  elastic  matrix)  developed  by  Rosen  (1965).  Davis  modeled  the  fibers  as  sinusoidal  in 
shape  and  found  shear  bucking  failure  to  be  responsible  for  significant  reductions  in  compression  strength. 
A  similar  model  was  developed  by  Hyer  (1986)  who  found  that  compression  strength  of  wavy  fibers 
decreased  significantly  with  decreasing  matrix  shear  modulus.  In  addition,  he  found  the  rumlinearity  of 
the  matrix  material  causes  a  further  reduction  in  compression  strength.  Fiber  kinking  failure  is  closely 
related  to  fiber  waviness.  Consequently,  numerous  investigations  have  been  conduaed  on  the  kinking 
failure  of  initially  misaligned  fibers  (Argon  1972;  Budiansky  1983;  Steif  1990a,  1990b;  Martinez  et  al. 
1981). 

The  influence  of  fiber  waviness  on  the  effective  elastic  mechanical  properties  of  composite  lamina  was 
investigated  by  Jortner  (1984).  The  paper  also  provides  a  thorough  review  of  the  literature  in  this 
area.  Jortner  discusses  the  typical  bounding  techniques  often  employed  in  predicting  the  influence  of  fiber 
wavir^ss  on  the  effective  elastic  properties  of  composites.  He  presents  an  exact  solution  for  predicting 
the  elastic  and  thermal  constants  of  wrinkled  regions  in  composites;  however,  the  model  is  limited  to 
waviness  that  exhibits  undulation  which  is  completely  in  phase  with  the  surrounding  material. 
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More  recently,  investigators  tuve  began  studying  the  effects  of  jdy  waviness  in  roultidireciionat 
laminates.  These  studies  have  focused  on  the  interactions  within  and  between  adjacent  [^es  within  a 
laminate  as  opposed  to  interactions  between  individual  fibers  and  the  matrix  material.  Shuait  (198S) 
studied  the  behavior  of  multidirectional  laminates  in  uniaxial  ctnnpression  emidoying  a  microbuckling 
model.  His  models  included  a  linear  buckling  analysis  for  ^it  waveloigth  buddng  and  a  geometrically 
nonlinear  analysis  for  wavy  layers.  He  observed  that  imeilaminar  shear  presses  were  a  maximum  su  a 
point  in  the  wavy  ply  where  die  out-of-(dane  undulation  is  greatest  0.e.,  the  point  of  inflection).  In 
addition,  he  concluded  that  laminate  failure  was  caused  by  either  short-wavelength  bucking,  imeilaminar 
shear,  or  in-fdane  shear  failures,  depending  on  ply  orientations  and  the  degree  of  waviness.  Total  laminate 
compressimi  strengdi  reductions  of  up  to  36%  were  reported.  Shuart  (1989)  subsequently  included  a 
matrix  onniMession  failure  mode  and  in-plane  ply  waviness  in  his  analysis.  His  work  was  further 
extended  by  Peel  (1990)  who  considered  additional  failure  modes  as  well  as  a  variety  of  wavy  ply 
geometries  and  configurations. 

Hyer,  Maas,  and  Fuchs  (1988)  and  Telegadas  and  Hyer  (1990)  have  studied  the  effects  of  ply  waviiKss 
(HI  the  stress  states  in  hydrostatically  loaded  cylinders.  They  employed  finite  element  analyses  to  study 
the  influence  the  geometry  of  the  waviness  and  material  properties  on  the  stress  states  that  develop.  They 
Observed  significant  interlaminar  stress  levels  in  wavy  plies  and  found  material  property  variations  to  be 
important  factors  in  the  magnitude  of  stresses  produced. 

The  influence  of  ply  waviness  m  the  performance  of  cylinders  has  also  been  investigated 
experimentally  by  Garala  (1989).  In  this  study,  numerous  8-in  morwcoque  graphite/epoxy  cylinders  were 
hydrostatically  loaded  to  failure.  Cylitaler  failures  might  have  been  initiated  by  interlaminar  shear  failure 
due  to  ply  waviness  in  die  hoop  and  axial  directions. 

The  approacdi  taken  in  diis  work  foUows  closely  that  of  Ishikawa  and  Chou  (1982,  1983a,  1983b, 
1983c)  in  their  studies  of  the  thermomechanical  behavior  of  two-dimensional  woven  fabric  composites. 
As  widi  the  models  they  developed,  the  theoretical  basis  for  the  model  developed  here  is  laminated  plate 
theory.  The  success  of  their  models  for  predicting  the  effective  thermomecdianical  properties  of  woven 
fabric  composites  jHOvides  initial  confidence  in  and  a  point  of  departure  for  the  results  generated  in  this 
work. 
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1.3  Problem  Statement  The  objective  of  this  woilc  is  to  quantify  the  influence  of  manufacturing 
induced  ply  waviness  on  the  stiflhess  and  strength  in  a  composite  and  to  assess  die  affect  on  stnictural 
performance.  Ply  waviness  is  assumed  to  exist  in  one  direction  only.  A  schematic  illustrating  ply 
waviness  is  shown  in  Figure  1,  where  the  wavy  ply  is  undulating  in  the  x-dircction.  This  layer  is  defined 
as  a  [0]  ply  which  is  embedded  in  [90]  plies  or  lamina  having  straight  fibers  aligned  in  the  y-direction. 

The  models  developed  in  this  woilc  piedia  the  in-plane  (x-y  plane)  longitudinal  mcxlulus, 
transverse  modulus,  Ey  in-plane  Poisson  ratio  and  shear  modulus,  and  respectively,  as  a  function 
of  the  undulation  geometry  and  material  properties  of  the  composite  lamina.  Tt^  in-plane  thermal 
expansion  coefficients,  and  (x^.  are  also  predicted. 

In  addition  to  the  elastic  properties,  the  model  developed  calculates  the  stresses  in  the  [0]  and  [90] 
plies  as  a  function  of  geometric  location  and  effective  in-plane  plate  loadings.  The  stresses  in  the  [90] 
ply  will  be  in-plane  only,  while  out-of-plane  stresses  in  the  [0]  {My  will  develop  due  to  the  inclination  of 
the  ply  out  of  the  x-y  plane.  The  magnitude  of  the  stress  components  in  each  ply  will  vary  along  the 
length  of  the  undulation  due  to  a  constant  in-plane  stress  resultant  assumption.  Once  the  stresses  in  the 
individual  plies  are  calculated,  the  maximum  stress  failure  criterion  is  applied  to  predict  ultimate  failure 
of  the  wavy  ply  configuration.  The  maximum  stress  failure  criterion  is  chosen  because  it  identifies  the 
failure  mode  and  provides  the  greatest  physical  insight  into  the  mechanisms  of  strength  reduction.  With 
the  failure  analysis  capability,  the  influence  of  waviness  on  strength  reduction  is  quantified.  Thermal, 
uniaxial,  and  biaxial  loadings  are  also  investigated. 

2.  ANALYSIS 

2.1  Background:  Laminated  Plate  Theory.  The  theoretical  development  of  the  analysis  is  presented 
in  this  section.  The  basis  for  the  analytic  model  developed  is  laminated  plate  theory  (Whitney,  Daniel, 
and  Pi|)es  1982;  Whitney  1987).  Laminated  plate  theory  is  well  develofied  and,  consequently,  only 
equations  pertinent  to  the  analysis  developed  here  are  presented. 

Laminate  res(X}nse  is  generally  described  through  a  relationship  between  effective  plate  loadings  and 
plate  deformation.  Typically,  the  relationship  between  in-plane  stress  and  moment  resultants,  N-  and 
A/,  ,  and  in-plane  strains  and  curvatures,  and  k,  .  is  given  by 
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(i.J  •  1,2,6) 


(1) 
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where  the  A^j,  Bfj,  and  D^j  stiffoess  matrices  are  defined  in  the  usual  maoier  as 


(Qij)^dz 


(2) 


or  simply 

Aij-  E  ('■»-*»-!) 
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8ii-E 
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The  lamina  stiffness  matrix  for  the  ply  is  iQip^,  n  is  the  total  number  of  plies  in  the  laminate,  and 
is  the  distance  from  the  laminate  midplane  to  the  ply  surface.  An  alternative  fonn  of  Equation  1  is 
obtained  through  its  inversion  and  is  given  as 
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where  Oij,  by,  and  dy  are  the  compliance  matrices  that  describe  the  effective  {riate-like  behavior  of  the 
laminate.  Ttrese  compliance  components  are  the  terms  ultimately  used  to  compute  the  effective  in-plaire 
thermomechanical  properties  of  the  laminate. 
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2.2  Description  of  the  Geometry.  In  order  to  quantify  the  influence  of  ply  waviness  tm  stiffiiess  and 
strength  reduction,  a  mathematical  description  of  the  geometry  is  required.  From  the  wavy  {^y 
configuratkm  depicted  in  Figure  1,  a  unit  cell  is  idealized  aiKl  is  shown  in  Figure  2.  A  half-sine  wave 
segmou  of  die  undulation  is  selected  as  a  representative  element  for  the  wavy  ply  configuratitm.  The  unit 
cell,  having  a  total  thickness  h,,  consists  of  an  unduling  [0]  ply  with  thickness  /flying  between  two  [90] 
regions.  The  amplitude  of  the  undulation  is  A.  The  half-wavelength  of  the  undulation  is  L  The  [90] 
plies  are  aligned  in  the  y-direction  and  are  assumed  to  have  no  undulation. 

Detailed  matiiematical  (tescriptions  of  the  individual  plies  within  the  unit  cell  expressed  as  a  function 
of  jc  are  required.  Referring  to  Figure  2,  the  following  distances  along  the  half-sine  wave  sluqied 
undulation  are  defined  by 


hi(x)  -  -h,n 

h^{x)  -  -All-hfll  f  [1  +  sin(l(x-i:))]  A/2 

fi^(x)  •  -AI2*hfl2  ♦  [1  *  sinilix-L))]  A/2 

f.  2 

h^{x)»h,l2  (5) 


A  local  out-of-plane  undulation  angle.  0,  is  also  needed.  Referring  to  Figure  3,  0  can  be  defined  as  a 
function  of  X  as 


e(x) 


■  tan~^ 

dih^ix)} 

*  tun 

dihjix)) 

dx 

dx 

(6) 


2.3  Effective  Elastic  Properties.  The  in-plane  thennomechanical  properties  of  tire  unit  <reil  are 
representative  of  the  re^nse  to  be  expected  of  the  entire  wavy  ply.  The  approach  taken  here  is  to  divide 
the  unit  cell  into  discrete  segments  of  width  dx  in  tire  direction  of  undulation  (see  Figure  3).  The  unit  cell 
is  defined  by  the  thickrress  of  tire  smallest  repeating  sublaminate  employed  during  fabrication,  as  well 
as  hf,  A,  and  L  defiired  previously.  A  systematic,  laminated  plate  theory  analogy  is  applied  to  eadt 
segment.  The  effective  response  of  the  unit  cell  is  obtained  through  an  averaging  approach  where  an 
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Figure  2. 


in-plane  stress  resultant  is  assumed  constant  widiin  each  discrete  segment.  In  this  a|^>roach.  average 
compliances  are  obtained  through  an  integration  procedure  and  subsequent  calculations  are  made  to 
evaluate  the  effective  thermomechanical  properties  of  the  unit  cell  as  a  function  of  ply  waviness. 


The  Aij,  Bfj,  and  D,y  stiffness  matrices  for  each  discrete  segment  of  width  dx  within  the  unit  cell  are 
first  computed  throu^  Equation  3  using  the  hfx)  values  defined  in  Equation  S.  The  unit  cell  is  essentially 
a  three-ply  laminate  with  x-dependent  A^j,  and  D^j  stiffness  matrices.  The  thickness  of  individual  piles 
are  defined  by  the  (90/0/90]  laminate  stacking  sequence  to  be  investigated.  The  Q^jS  for  the  [90]  plies 
are  constant  along  the  x-direction  and  are  computed  from  the  typical  engineering  constants  of  tire  straight 
lamina  (i.e.,  £j,  £2.  Vjj.  and  ^12)  in  the  usual  fashion  (Whitney,  Daniel,  and  Pipes  1982;  Whitney  1987). 
The  Q^jS  for  the  (0]  ply,  however,  must  take  into  account  the  out-of-plane  (z-direction)  orientation.  The 
effective  engineering  properties  used  to  compute  the  Q^s  in  each  dx  segment  along  the  x-direction  are 
a  function  of  the  inclination  angle,  6,  and  are  given  by 


(®)  -  +  (1/G,3  -  2vi3/£i)/e^Wo^  + 

Ey(6)  -£2 


Vx,(0) 

GxyO) 
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^e''l2  ^  ''32'”0 
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Gi2 

(q  »  cosd 


mg  »  Sind 


(7) 


where  £j,  £3,  Gj2.  V12,  and  G23  are  the  typical  principal  system  engineering  constants  of  the  straight 
lamina.  The  effective  a,y,  b,y,  and  d,y  compliance  matrices  for  each  dx  segment  within  the  unit  cell  are 
computed  through  inversion  of  the  Aij,  Bij,  and  D^j  stiffness  matrices  according  to  Equation  4. 

The  local  warping  constraint  is  imposed  on  each  dx  segment  in  the  unit  cell.  This  assumption  leads 
to  a  more  realistic  estimate  of  the  effective  response  of  the  wavy  ply  configuration  since  no  global 
out-of-plane  warping  is  expected  in  an  axisymmetric  geometry  or  in  a  laminate  of  appreciable  thickness. 
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Under  the  local  warping  constraint,  a  zero  curvature  and  non-zero  moment  condition  are  imposed.  The 
bending  free  or  local  warping  constraint  behavior  of  each  dx  segment  can  be  expressed  as  a  function  of 
the  full  a;j,  bij,  and  d^j  compliance  matrices  defined  in  Equation  4  (Ishikawa  and  C3iou  1983a), 

a'j  .  (8) 


where  a,y*  is  now  the  in-plane  compliance  matrix  with  local  waiping  constrained. 

In  order  to  calculate  the  average  or  effective  response  of  the  wavy  ply,  an  average  in-plane  compliance 
matrix  is  defined.  Assuming  a  uniform  in-plane  stress  resultant  exists  in  eadi  dx  segment,  the  average 
in-plane  compliance  matrix  is  obtained  through  the  following  integration  over  the  entire  unit  cell  lengdi: 


Here  a,y**  is  the  effective  in-plane  compliance  matrix  for  the  wavy  ply  with  the  local  waiping  condition 
imposed.  Similarly,  the  compliance  matrices  b-j  and  d^j  could  be  defined,  however,  these  matrices  are  not 
employed  in  determining  the  effective  in-plane  mechanical  properties  of  the  wavy  fiy  construction  (see 
Equation  8). 


The  effective  in-plane  mechanical  properties  of  the  wavy  ply  are  computed  from  the  following  a-** 
compliance  coefficients  according  to 
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It  is  noted  that  the  effective  in-plane  mechanical  propeities  with  local  waiping  allowed  could  have 
beat  computed  in  a  more  straightforward  fashion  than  those  computed  with  the  local  waiping  constrained 
condidcKi  applied.  The  procedure  is  essentially  identical  to  that  outlined  above  with  the  step  taken  in 
Equation  8  omitted.  This  procedure  would  predict  stiffnesses  significantly  more  compliant  due  to 
unconstrained  waiping  of  the  unit  cell. 

2.4  Effective  Thermal  Expansion  Coefficients.  A  similar  approach  for  the  prediction  of  the  effective 
mechanical  properties  of  the  wavy  ply  construction  is  taken  in  calculating  the  effective  thermal  expansion 
coefficients.  In-plane  thermal  load  resultants.  N^,  in  each  dx  segment  of  the  unit  cell  are  required  and 
are  computed  through  the  following  equation: 

where  AT  is  the  unifonn  theimal  load,  and  the  ct-s  are  the  in-plane  thennal  expansion  coefficients  of  the 
[0]  and  [90]  plies.  Here  again  both  and  ay  in  the  [90]  ply  are  constant  along  tlK  x-direction  within 
the  unit  cell.  The  [0]  ply  properties  are  functions  of  ^(x)  (Equation  6)  and  are  transformed  into  the 
x-y  coordinate  system  from  the  principal  coordinate  system  defined  by  the  local  fiber  orientation  shown 
in  Figure  3.  Thermal  expansion  coefficients  for  each  segment  are  defined  by 

where  is  that  computed  in  Equation  8.  The  effective  in-plane  thennal  expansion  coefficients  for  the 
wavy  ply  configuration,  a*,*  are  then  calculated  by  the  same  averaging  approach  described  above  through 

This  completes  the  theoretical  development  of  the  analytic  model  for  predicting  the  effective 
thermomechanical  propeities  of  the  wavy  ply  configuration.  In  the  following  section,  the  determination 
of  the  individual  [0]  and  [90]  ply  stresses  and  strength  analysis  are  discussed. 

2.5  Ply  Stresses.  In  order  to  gain  a  fundamental  understanding  of  the  failure  mechanisms  associated 
with  wavy  {dies  and  to  quantify  the  influence  of  the  degree  of  undulation  on  strength  reduction,  individual 
ply  stresses  in  the  wavy  ply  configuration  must  be  computed.  Individual  ply  stresses  in  both  the  [0]  and 
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[90]  plies  are  computed  from  given  or  defined  in-plane  stress  resultants  in  the  jc-  and  y-dircctions 
(Figure  1).  Recalling  that  a  constant  in-planc  stress  resultant  within  the  unit  cell  was  assumed,  the 
individual  ply  stresses  are  expected  to  vary  along  Lhe  x-directlon  within  the  unit  cell  due  to  the  undulating 
[0]  ply.  The  procedure  for  computing  these  stresses  is  outlined  below. 

For  a  given  set  of  in-plane  stress  resultants  due  to  thermal  or  mechanical  loading,  Nj,  the  in-plane 
strains  within  each  dx  segment  within  the  cell  are  computed  by 

k}-k]h) 

where  a*  is  the  compliance  matrix  computed  in  Equation  8.  It  should  be  understood  that  the  in-plai^ 
strains  (defined  in  the  x-  and  y-directions)  of  each  ply  within  a  single  dx  segment  are  identical  throu^ 
the  thickness  (z-direciion)  but  vary  along  the  length  (x-direcdon).  Since  local  waiping  is  constrained,  tl» 
in-plane  strains  represent  the  total  ply  strains  and  the  typical  strain  displacement  relation  is  not  required. 
The  individual  ply  stresses  (in  the  x-  and  y-directions)  within  each  dx  segment  are  evaluated  through 
Hooke’s  law  relation  given  by 

{(o’)*}  -[(fiy)*  ]{(e;)*-(ap*Ar} 

where  the  Q^jS  are  the  same  as  those  used  in  the  effective  mechanical  property  calculations. 

While  Equation  IS  provides  the  ply  stresses  in  each  dx  segment  defined  in  the  global  x-  and 
y-directions,  these  stresses  must  be  transfonned  into  the  principal  system  of  the  [0]  and  [90]  plies.  The 
principal  material  coordinate  system  for  a  unidirectional  lamina  is  shown  in  Figure  4.  Referring  to 
Figures  1  and  4,  the  longitudinal  or  1 -direction  stresses  in  the  [90]  plies  will  correspond  to  y-direction 
stresses  calculated  in  Equation  15.  Transverse  or  2-direction  stresses  in  the  [90]  plies  will  correspond  to 
the  calculated  x-direction  stresses. 

Referring  to  Rgures  1  and  4,  the  2-direction  stresses  in  each  dx  segment  of  the  [0]  ply  correspond  to 
the  calculated  y-direction  stresses.  The  situation  is  slightly  more  complicated,  however,  for  computing 
the  other  principal  stress  components  in  the  [0]  ply.  Since  the  [0]  ply  is  undulatii^  out  of  the  x-y  plane, 
transverse  normal  (3-direction)  and  interlaminar  shear  stress  (1-3)  components  will  develop  within  the  unit 
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Figure  4.  Lamina  Cooidinate  System. 

celL  Theae  stress  components  result  solely  frcnn  the  x-direction  stresses  and  are  computed  through 
tensotal  transfonnation  of  the  x-diiection  stress  components. 

In  summary,  the  [90]  {riles  will  contain  1*  and  2-direction  stress  components  only.  The  [0]  {riy  will 
coittain  1,  2,  3  and  1-3  stress  com(X)nents  due  to  the  out-of-plane  undulation  within  the  unit  well. 
Furthermore,  these  stress  compxments  will  vary  along  the  x-direction  within  the  unit  cell.  Understanding 
these  stress  variations  will  (xovide  fundamental  insight  into  the  dominant  failure  mechanisms  characteristic 
of  wavy  {riies.  In  the  next  sectkHi,  die  failure  criterion  ap{riied  to  the  wavy  ply  configuraticm  is  firesented. 

2.6  Applied  Failure  Criterion:  Maximum  Stress.  The  failure  analysis  developed  in  this  investigation 
is  iMsed  on  the  maximum  stress  failure  criterion.  Failure  of  the  wavy  {riy  configuration  is  defined  when 
a  prindpal  stress  omponem  in  either  the  [0]  or  [90]  {riies  of  any  dx  segment  altmg  die  unit  ceU  reaches 
its  corresponding  strength  allowrririe.  To  provide  results  that  are  useful  in  making  design  decisions  based 
on  the  failure  analysis,  a  scaling  a{)proach  which  determines  the  effective  in-plane  load  resultants  at  failure 
was  adt^ried.  For  a  given  wavy  ply  geometry  configuration,  the  in-plane  load  resultants  for  each  dx 
element  along  the  unit  ceil  which  would  initiate  failure  are  computed.  The  lowest  load  rnultartt  for  any 
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segment  over  the  entire  length  of  the  unit  cell  defines  the  loads  at  failure  for  the  entire  wavy  ply 
configuration.  Strength  reduction  associated  widi  ply  waviness  is  based  on  comparisons  between  failure 
loads  of  the  wavy  ply  configuration  to  those  of  a  similar  straight  laminate  configuration  having  zero 
amplitude  of  undulation. 

Bodi  in-{dane  medianical  loading  and  thermal  loading  have  been  incorpontted  into  the  analysis. 
Thermal  loading  is  applied  in  the  form  of  a  uniform  temperature  differential.  Themud  loads  at  ftulure  ate 
scaled  in  the  same  manner  described  above.  In  addition,  a  combination  of  thermal  and  mechanicad  loading 
as  well  as  biaxial  loading  ate  permitted. 

3.  RESULTS 


Parametric  studies  illustrating  the  results  from  the  analytical  models  developed  above  for  predicting 
the  thermomechanical  properties,  ply  stresses,  and  strength  reduction  in  wavy  plies  are  presented  in  this 
section.  Significant  findings  are  repotted  in  the  body  of  this  report  while  less  interesting  results  generated 
have  been  documented  in  the  Appendix.  Before  the  results  are  presented,  a  constraim  tm  the  geometric 
parameters  describing  the  wavy  ply  configuration  is  discussed.  Results  are  generated  for  both 
AS4-Carbon/PEKK  and  S2-Glass/PEKK  composite  material  systems.  A  summary  of  the  ge:>metric 
parameters,  material  properties,  and  strength  allowable  used  in  the  parameuic  studies  is  presented. 


3.1  Geometric  Constraint  on  Plv  Waviness.  To  describe  the  physical  geometry  of  the  wavy  ply 
configuration,  a  parameter  that  identifies  the  ratio  of  the  total  [90]  ply  thickness  (i^  to  [0]  ply  thickness 
(to),  TX  (thickness  ratio),  is  required 


2  'tgQ  ^  2  "tgQ 


(16) 


The  definition  of  this  ratio  is  illustrated  in  Figure  5.  Once  a  TX  ratio  is  defined,  the  total  thickness 
of  the  unit  cell  becon^ 

A,  -Ay(l+rR)  (17) 


A  geometric  constraint  imposed  on  the  amplitude  of  the  [0]  ply  undulation  wavy  ply  configuration  for  a 
given  A,  and  is  identified  in  Figure  6  to  be 


A  <  hf  -  hf 


(18) 


In  the  remainder  of  this  section,  parametric  studies  are  presented  to  illustrate  the  ngnificant  findings 
of  this  woik.  Presentation  of  the  results  is  facilitated  with  the  comi^ehensive  iOustration  indicating  all 
details  of  the  wavy  ply  configuration  including  the  global  and  ply  coordinate  systems  used  in  the  analysis 
(see  Figure  7). 

3.2  Parametric  Studies.  Parametric  studies  were  conducted  between  zero  anqiiitude  and  the  value 
defined  in  Equation  18  to  bound  the  potential  process-induced  reduction  in  stiffiiess  and  strength. 
Photomicrographs  of  actual  filament  wound  composite  cylinders,  however,  revealed  that  the  nominal 
unplitudes  were  much  less  than  this  ut^r  bound.  Photomicrographs  were  useful  in  defining  a  range  of 
half-wavelengths  to  be  investigated  in  this  study.  The  most  ounmon  wavelength  observed  was 
afqaoximately  20  ply  thicknesses,  which  correlated  well  with  the  nomiiud  width  of  the  tow  employed 
during  filament  winding.  The  geometric  parameters  studied  are  listed  in  T^e  1.  A  constant  [0]  {Hy 
thickness  of  0.(X)S  inches  was  used  as  the  baseline.  Tables  2  and  3  summarize  the  ejqrerimentally 
determined  mechanical  properties  and  strength  allowables  of  the  AS4  CarborVPEKK  and  S2  Glass/PEKK 
composite  systems  used  in  the  study. 

3.2.1  Effective  Elastic  Properties.  The  influence  of  various  geometric  parameters  on  the  effective 
elastic  properties  of  the  wavy  ply  configuration  were  investigated  first  The  influence  of  undulation 
amplitude,  A,  on  the  reduction  of  j;-direction  modulus,  £^,  for  the  AS4  Graphiie/PEKK  composite  is  shown 
in  Figure  8  for  various  TR  ratios  and  a  unit  cell  lengdi  of  L  =  5*)y.  Stiffness  rcducticm  is  defined  as  the 
ratio  of  wavy  ply  stiffness  to  the  stiffness  of  an  equivalent  laminate  construction  widi  a  zero  undulation 
amplitude.  The  results  indicate  significant  stiffness  reduction  of  up  to  60%  for  A//iy  values  below  5  in  this 
extreme  case.  The  influence  of  increasing  the  unit  cell  length  to  £  =  20*/y  is  shown  in  Figure  9.  Stiffness 
reduction  was  not  as  severe  (<30%)  for  A//y  values  below  5.  The  same  study  for  a  unit  cell  lengdi  of 
L  =  5»hf  for  the  S2  Glass/PEKK  composite  is  shown  in  Figure  10.  This  study  indicates  that  stiffness 
reductions  in  the  glass  system  are  not  as  sensitive  to  undulation  amplitude  as  in  the  cartwn  system.  This 
is  true  because  the  graphite  system  exhibits  substantially  more  anisotro|diy  in  the  fiber  direction  than  the 
glass  composite.  Other  in-plane  mechanical  properties  such  as  transverse  stiflhess,  Ey,  poisson  ratio, 
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?  [90]  Ply 


[0]  Ply 


Figure  7.  Global  and  Plv  Coordinate  Systems  Used  in  Analysis. 
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(0=V)  *3  /  *3 
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(0=V)  *3  /  *3 


Figure  10.  InfluCTce  of  Undulation  Amplitude  cti  in  S2  Glass/PEKK  5). 


20 


Table  1.  Parameter  Study:  Geometry 


Table  2.  Material  Properties 


Property 

AS4  Graphite/PEKK 

S2  Glass/PEKK 

(Msi) 

17.0 

8.0 

E2  (Msi) 

1.40 

2.75 

Gj2  (Msi) 

0.98 

1.01 

G23  (Msi) 

0.98 

0.59 

''12 

0.36 

0.27 

''23 

0.36 

0.31 

aj  (pe/®F) 

-0.5 

2.3 

02(Me/“F) 

15.5 

18.5 

a3(p£/"F) 

15.5 

18.5 

Table  3.  Strength  Allowables 


Allowable 

AS4  Graphite/PEKK 

S2  Glass/PEKK 

Xir(ksi) 

28S 

243 

X2r(ksi) 

5.2 

XSrOcsi) 

5.2 

8.5 

XIC  (ksi) 

155 

177 

X2C(ksi) 

26.5 

30.6 

X3C  (ksi) 

26,5 

35.0 

S13  (ksi) 

15.2 

17.0 

S12  (ksi) 

19.0 

15.7 

in-plaire  shear  modulus,  <7^.  were  found  to  be  relatively  insensitive  to  undulation  amplitude  in  the  x- 
direction  for  both  material  systems  studied  (see  Appendix  A). 

In  general,  stiffness  reduction  is  most  severe  in  the  carbon  system,  increasing  with  increasing 
undulidon  amplitude  and  decreasing  unit  cell  length. 

3.2.2  Effective  Thermal  Expansion  Coefficients.  Increasing  undulation  amfditude  was  found  to 
increase  the  ^'direction  thermal  expansion  coefficient  up  to  70%  in  the  carbon  system  with  L  =  20</y 
(figure  11).  A  30%  increase  in  the  glass  system  is  predicted  for  L  =  (figure  12).  The  increase  is 
associated  with  the  out-of-plane  rotation  of  the  axial  plies  which  decreases  axial  stiffness  and  provides  less 
constraint  on  thermal  deformation  of  the  [90]  plies.  Changes  in  transverse  thermal  expanson  were  found 
to  be  insignificant  in  all  cases  investigated  (see  Appendix  A). 

3.3  Ply  Stresses.  Ply  stresses  along  the  unit  cell  length  predicted  from  the  analysis  are  presented  to 
provide  insight  into  the  dominant  failure  mechanisms  characteristic  of  the  wavy  |Hy  configuration.  Ply 
stresses  are  nondimensionalized  by  the  corresponding  strength  allowables  presented  for  the  AS4 
Graphite/PEKK  composite  system  in  Table  3.  The  magnitude  of  the  ply  stresses  are  calculated  from  a 
uniaxial  x-direcUon  load  resultant  that  would  cause  failure  of  an  identical  laminate,  except  with  a  zero 
amplitude  wavy  ply  configuration  -  901  IbAm).  Ply  stresses  in  the  [0]  layer  in  the  principal  1 -direction 
are  presented  in  figure  13  for  undulation  amplitudes  A  <,  0.01  in.  Stress  components  are  shown  to  be  a 
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Figure  11.  Influence  of  Undulation  Amplitude  cm  in  AS4  Graphite/PEKK  (L/Ay«  20). 
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Figure  12.  Influence  of  Undulation  Amplitude  on  in  S2  Glass/PEKK  (L/A^=  5). 
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minimum  in  regions  within  the  unit  cell  where  out-of-plane  orientations  are  maximum.  In  this  situation, 
the  load  is  transferred  to  the  2-direction  of  tte  [90]  plies  to  satisfy  force  equilibrium. 

Interlaminar  normal  (3-direction)  and  shear  stress  components  (1-3)  in  the  [0]  jdy  ate  shown  in 
Figures  14  and  15  for  the  same  load  situation.  Normal  stress  components  are  ngniflcam  in  the  region  of 
maximum  undulation,  teaching  just  beyond  the  strength  allowable  in  this  example.  The  interlaminar  shear 
stress  componoits,  however,  are  far  more  significant  teaching  nearly  three  limes  the  strength  allowable 
for  this  extreme  case.  Other  less  significant  results,  including  [90]  ply  stress  {dots  and  load  cases  in  the 
y-diiection,  are  documented  in  Appendix  B. 

The  most  significant  finding  of  the  parametric  study  is  the  high  interlaminar  shear  stress  compoKnts 
that  develop  in  the  [0]  plies  due  to  the  out-of-plane  orientation.  The  stress  levels  observed  indicate  that 
this  is  a  dominant  failure  mode  when  the  wavy  ply  configuration  is  loaded  in  the  direction  cf  undulation. 
It  also  quantifies  the  influence  of  undulation  amplitude  on  the  ply  stress  levels  that  develop  and 
demonstrates  the  importance  of  controlling  ply  waviness  in  the  load  direction  during  manufacture.  These 
observations  are  also  consistent  with  those  made  by  Hyer,  Maas,  and  Fuchs  (1988)  and  Telegadas  and 
Hyer  (1990). 

3.4  Strength  Reduction  Due  to  Plv  Waviness. 

3.4.1  Uniaxial  Loading.  Gearer  insight  into  the  influence  of  ply  waviness  <m  strength  reduction  is 
achieved  by  examining  the  reduction  in  the  ultimate  applied  load  resultant  that  the  wavy  ply  configuration 
can  withstand  before  the  stress  component  in  any  ply  reaches  its  corresponding  strength  allowable.  In 
Figures  16  and  17,  strength  reduction  in  AS4  Graphite/PEKK  wavy  ply  configurations  of  (A  ^.01  in) 
length  L  =  5*hf  and  20*/y  are  presented  for  a  uniaxial  loading  in  the  x-direction.  These  results  indicate 
that  significant  strength  reduction  (65%  for  L  =  5»h^and  35%  for  L  =  20*/t^  along  the  unit  cell  length  is 
possible.  It  is  pointed  out  that  the  dominant  failure  mode  transitions  fmn  a  longitudinal  compressive 
failure  (XIQ  in  the  [0]  ply  to  interlaminar  shear  failure  (513)  in  the  [0]  ply  in  regions  where  the  out-of¬ 
plane  orientation  increases.  Other  results  ^ow  that  ply  waviness  in  the  x-direction  does  not  significantly 
affect  the  y-direction  strength  allowable. 
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gure  15.  Interlaminar  Shear  Stress  Within  a  Wavv  fOI  Plv  in  AS4  Graphite/PEKt’ 
Loading  (JR  =  l,Uh  =  5.  (/I  =  0)  =  -901  Ib/in  at  Failure). 


Uniaxial 


28 


x/L 


Figure  16.  Strength  Reduction  in  AS4  Graphite/PEKK:  Uniaxial  Loading  iTR  =  2.  Uh^  =  S. 
M  =  01  =  -901  Ib/in  at  Failure). 
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Figure  17.  Strength  Reduction  in  AS4  Graphite/PEKK:  Uniaxial  Lxiading  (TR  =  2.  L/hr-  20. 
M  =  0)  g  -901  Ib/in  at  Failure!. 
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3.4.2  Biaxial  Loading.  Nondimensional  failure  envelopes  for  biaxial  load  conditions  with  and  without 
thermal  loading  were  generated  to  quantify  the  influence  of  ply  wavincss  on  strength  reduction.  The 
uniaxial  stress  resultants  used  to  nondimensionalize  the  biaxial  failure  envelopes  presented  here  are 
(/i  =  0)  =  -901  Ib/in  at  failure  and  (A  =  0)  =  -1,612  Ib/in.  The  influence  of  ply  undulation  on  strength 
reducticMi  without  thermal  loading  is  illustrated  in  Figure  18.  The  plot  indicates  the  y-direction  strength 
is  insensitive  to  ply  waviness  while  x*direcdon  strength  can  be  substantially  reduced  (up  to  3S%  for  the 
particular  ply  configuration  investigated).  The  same  failure  envelope,  excein  with  a  thermal  load  of 
AT  =  -150®  F,  is  presented  in  Figure  19.  Since  the  thermal  loading  primarily  influences  the  transverse 
stresses  in  the  [90]  plies,  the  addition  of  the  thermal  load  has  little  affect  on  the  ultimate  strength  of  the 
applied  loads. 

The  ultimate  y-direction  loading  is  defined  almost  exclusively  by  the  longitudinal  compressive  strength 
allowable  (XIC)  of  the  [90]  plies.  In  contrast,  and  consistent  with  the  ply  stress  results  presented  earlier, 
the  x-direction  loading  results  in  interlaminar  shear  failure  (513)  in  the  [0]  ply.  This  failure  mechanism 
defines  the  right-hand  side  of  the  failure  envelopes  and  is  highly  sensitive  to  the  waviness  of  the  [0]  ply. 

A  major  conclusion  drawn  from  these  results  is  that  strength  reduction  in  the  direction  of  ply  waviness 
is  significant,  highly  sensitive  to  the  degree  of  waviness,  and  is  dominated  by  interlaminar  shear  failure 
in  the  wavy  ply.  Strength  reduction  in  the  non-wavy  direction  is  relatively  insignificant. 

4.  CONCLUSIONS 

An  analytical  model  was  developed  to  quantify  the  influence  of  ply  waviness  on  stiffness  and  strength 
reduction.  The  model  provides  an  efficient  design  tool  for  accessing  the  synergistic  effects  of  ply 
waviness  on  laminate  stiffness  and  strength  which  can  be  directly  related  to  structural  performance. 
Parametric  studies  were  conducted  to  identify  the  criticality  of  ply  waviness  on  the  stiffness  and  strength 
reduction  of  realistic  composite  structures. 

Several  significant  findings  were  made.  Stiffness  reduction  due  to  ply  waviness  is  most  significant 
in  the  direction  of  ply  undulation.  Other  in-plane  mechanical  properties  are  relatively  insensitive  to  ply 
waviness.  The  magnitude  of  stiffness  reduction  is  extremely  sensitive  to  the  degree  of  ply  waviness. 
Increasing  undulation  amplitude  and  decreasing  unit  cell  length  boih  increase  observed  stiffness  reduction. 
Greater  stiffness  reduction  can  be  expected  in  systems  with  higher  material  anisotropy  (i.e.,  «  1). 
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Figure  18.  Nondimensional  Biaxial  Compression  Failufe  Envelopes  in  AS4  Graphite/PEKK 

(TR  =  2.  Uhj  =  20.  ^T  =  0°  F.  N.  M  =  Ot  =  -901  Ib/in.  N^.  (A  =  Q)  =  -1.612  Ib/in). 
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Dominant  failure  mechanisms  governing  ultimate  failure  in  wavy  [dies  and  the  influence  of  ply 
waviness  on  strength  reduction  were  also  quantified.  The  most  significant  strength  reduction  results  from 
interiaminar  shear  failure  of  the  wavy  ply  when  loaded  in  the  direction  of  undulatioa  Strength  reduction 
is  extremely  sensitive  to  the  degree  of  ply  waviness.  Failures  associated  with  loading  transverse  to  the 
direction  of  undulation  were  unaffected  by  the  degree  of  ply  waviness.  Failures  associated  with  losing 
transverse  to  tire  directimi  of  undulation  were  unaffected  by  the  degree  of  ply  wavtness  and  were 
essentially  the  same  as  straight  laminate  failure  loads. 

The  analysis  developed  in  this  work  plays  an  impoitant  role  in  the  design  of  composite  structures 
where  ply  waviiress  cannot  be  ignored.  The  design  procedure  should  employ  the  predicted  reduced 
stiffness  properties  associated  with  the  manufactured  ply  waviness  into  a  suitable  global  stress  analysis. 
The  resulting  global  stress  predictions  should  then  be  incorporated  back  into  the  failure  analysis  developed 
in  this  work  to  provide  realistic  estimates  of  the  strength  reduction  to  be  expected.  Predicted  stiffiress  and 
strength  reduction  due  to  ply  waviness  represents  one  mechanism  that  may  explain  why  design  and  failure 
analysis  of  thick-seaicHied  composite  structures  based  on  flat,  autoclave-cured  coupon  property  data  do 
not  correlate  well  with  test  results. 
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1  David  Taytor  Research  Center 
Ship  Structures  and  Protection  DepL 
ATTN:  Dr.  J.  Corrado,  Code  1702 
Bethesda,  MD  20084 

4  DirecbM’ 

Lawrence  Livermore  Natkxial  Laboratcny 
ATTN:  R.  Christensen 
S.  deTeresa 
W.  Feng 

F.  Magness 
P.O.  Box  808 
Livermore,  CA  94550 

2  Pacific  Northwest  Laboratory 

A  Div  of  Battelle  Memorial  InsL 
Technical  Information  Section 
ATTN:  M,  Smith 
M.  Garnich 
P.O.  Box  999 
Richland.  WA  99352 

6  Director 

Sandia  National  Laboratories 
Applied  Mechanics  Dept., 

Division-8241 
ATTN:  C.  Robinson 

G.  Benedetti 
W.  Kawahara 
K.  rerano 

D.  Dawstm 
P.  Niclan 
P.O.  Box  969 

Livemtore,  CA  94550-0096 
1  Director 

Los  Alamos  National  Laboratory 
ATTN:  D.  R^m 
MEE-13,  Mail  Slop  J-576 
P.O.  Box  1633 
Los  Alamos,  NM  87545 

3  University  of  Delaware 
Center  for  Composite  Materials 
ATTN:  J.  Gillespe 

B.  Pipes 
M.  Santare 

201  Spencer  Laboratory 
Newark.  DE  19716 
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North  Carolina  State  University 
Civil  Engineering  Dept. 

ATTN:  W.  Rasdorf 
L.  Spainhour 
P.O.  Box  7908 
Raleigh,  NC  27696-7908 

University  of  Utah 
Department  of  Mechanical  and 
Industrial  Engineering 
ATTN:  S.  Swanson 
Salt  Lake  City,  UT  84112 

Stanford  University 
Department  of  Aeronautics  and 
Aeroballistics  Durant  Bldg. 
ATTN:  S.  Tsai 
Stanford,  CA  94305 

Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 
ATTN:  Richard  McNitt 
227  Hammond  Building 
University  Park,  PA  16802 

Pennsylvania  State  University 
ATTN:  Rena  S.  Engel 
245  Hammond  Building 
University  Park,  PA  16802 

Pennsylvania  State  University 
ATTN:  David  W.  Jensen 
233-N  Hammond 
University  Park,  PA  16802 

Alliant  Techsystems,  Inc. 

ATTN:  C.  Candland 

J.  Bode 

K.  Ward 
5640  Smetana  Dr. 

Minnetonka,  MN  55343 


1  University  of  Illinois  at 

Urbana-Champaign 
National  Center  for  Composite 
Materials  Research 
216  Talbot  Laboratory 
ATTN:  J.  Economy 
104  South  Wright  Street 
Urbana.IL  61801 

2  Olin  Corporation 
Flinchbaugh  Division 
ATTN:  E.  Steiner 

B.  Stewart 
P.O.  Box  127 
Red  Lion,  PA  17356 

1  Olin  Corporation 

ATTN:  L.  Whitmore 
10101  9th  Sl.  NotUi 
Sl  Petosburg,  FL  33702 

1  Chamberlain  Manufacturing  Corporation 

Research  and  Development  Division 
ATTN:  T.  Lynch 
550  Esther  Street 
P.O.  Box  2335 
Waterloo,  lA  50704 

1  Custom  Analytical  Engineering  Systems,  Inc. 

ATTN:  A.  Alexander 

Star  Route,  Box  4A 
Rintstone,  MD  21530 

2  Institute  for  AdvatKed  Technology 

ATTN:  T.  Kiehne 

H.  Fair 

4030-2  W.  Braker  Lane 
Austin,  TX  78759 

1  UCLA 

MANE  Department,  Engineering  IV 

ATTN:  H.  T.  Hahn 

Los  Angeles,  CA  90024-1597 

2  Virginia  Polytechnical  Institute  and 

State  University 
DepL  of  ESM 
ATTN:  M.  W.  Hyer 

K.  L.  Reifsnider 
Blacksburg,  VA  24061-0219 
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2  University  of  Dayton  Research  Institute 
ATTN:  R.  Y.  Kim 
A.  K.  Roy 

300  Ct^lege  Park  Ave. 

Dayton.  OH  4S469-0168 

1  Universily  of  Dayton 

ATTN:  J.M.  Whitney 
300  Ctdiege  Park  Ave. 

DaytMi.OH  45469-0240 

1  University  of  Kentucky 

ATTN:  L.  Penn 
763  Anderson  Hall 
Lexington.  KY  40S064X)46 

1  Purdue  University 

School  of  Aeronautics  and 
Astronautics  ,  ' 

ATTN:  C.  T.  Sun  ,  - 
West  Lafayei»,'IN'*47907-1282 

1  Drexetl  University 
ATTN:  Albert  S.D.  Wang 
32nd  and  Chestnut  Sl 
Philadelphia,  PA  19104 

2  Kaman  Sciences  Corporation 
ATTN:  D.  Elder 

J.  Betz 
P.O.  Box  7463 

Colorado  Springs.  CO  80933 


iNTCNTKmALLY  LEFT  BLANK. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  contiiuiing  effort  to  improve  the  quality  of  the  reports  it  puUishes.  Your 
conunents^^wers  to  the  itemaApestions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  AltL-TR-IlO _ DareofReport  April  1993 _ 

2.  Date  Report  Received _ 

3.  I>ces  this  report  satisfy  a  need?  (Cbmmem  on  purpose,  related  project,  or  other  area  of  iruerest  for 

which  the  report  will  be  used.) _ 


S.  Has  the  information  in  this  report  led  to  any  quantitruive  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Cotieoion,  please  provide  the  Current  or  Comxt  address 
above  and  the  Old  or  biconect  address  below. 

Organization 

OLD  Name 

ADDRESS _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 

(DO  NOT  STAPLE) 


Department  of  the  Army 


OFROAL  BUSINESS 


BUSINESS  REPLY  M.UIL 
RUST  BjBS  rew  It  aooi.  IW,  w 

Po3iagc  will  &«  paid  by  addiassM 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Ub) 
Aberdeen  Proving  Ground,  MD  21005-5066 


